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The Landau Fermi Liquid parameters are calculated for charge neutral asymmetric nuclear 
matter in beta equilibrium at zero temperature in the presence of a very strong magnetic 
field with relativistic mean-field models. Due to the isospin structure of the system, with 
different populations of protons and neutrons and spin alignment to the field, we find non- 
vanishing Landau mixing parameters. The existence of quantized Landau levels for the 
charged sector has some impact on the Landau parameters with the presence of discretized 
features in those involving the proton sector. Using the Fermi liquid formalism singlet and 
triplet excited quasiparticle states are analyzed, and we find that in-medium effects and 
magnetic fields are competing, however, the former are more important in the interaction 
energy range considered. It is found that for magnetic field strengths Logio B (G) < 17 
the relative low polarization of the system produces mild changes in the generalized Landau 
parameters with respect to the unmagnetized case, while for larger strengths there is a 
resolution of the degeneracy of the interaction energies of quasiparticles in the system. 
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I. INTRODUCTION 



Asymmetric nuclear matter is presently an important topic proposed for study in experiments at 
radioactive beam facilities such as FAIR [1] at GSI, SPIRAL2 [2] at GANIL, ISAC-III at TRIUMF 
[3] or FRIB [4] at MSU, among others [5]. These will allow the investigation of nuclei in regions of 
the nuclear chart far from the stability line. These nuclear regions, where the isospin asymmetry 
ratio, given by the ratio of the proton vector number density, rip, to baryonic vector number 
density, ng, defined as Yp = rip/riB, largely departs from the 0.5 value, are of interest in the study 
of stability of exotic nuclei as in the neutron rich nuclei. Recently studied examples of those are the 
isotopes of Ca [6] with Z = 20 or the isotones N = 30 of Ca and Sc with CLARA experiment [7] in 
Legnaro. The study of observables related to nuclei far from the isospin stability line allows for an 
improved description of neutron rich environments as those of interest in astrophysical scenarios 
of neutron star matter. They are relevant to interiors of compact stellar objects like neutron stars 
(NS) arising in the aftermath of a supernova event. In this context, the rapid deleptonization of 
the NS following the gravitational collapse of the inner regions with the proton-electron capture 
process makes matter more neutron rich since (anti) neutrinos diffuse out of the star [8, 9]. Also 
supernova matter can be considered to be constitued by a set of nuclei, in the nuclear statistical 
equilibrium (NSE) approximation, and presents a distribution of nuclei shifted to the N > Z region 
[10, 11]. 

On Earth, the high temperature (T) and small baryonic chemical potential (^_b) region of 
matter phase space has been somewhat tested [12] at RHIC, and it will be possible to further 
test with other heavy-ion experiments like ALICE [13] at CERN. Although this improvement of 
our knowledge of the phases of matter is certainly valuable, the phase diagram of nuclear matter 
relevant for the equation of state (EOS) of NS is that of cold asymmetric nuclear matter in the low 
temperature range. 

Historically, most of the existing literature on the nuclear matter EOS characterization [14] 
has neglected the input from external fields, in particular, in the case of magnetic fields, this is 
claimed due to the tiny value of the nuclear magnetic moment [15]. On terrestrial experiments 
recent estimations of the magnetic field strength that could be produced dynamically at CERN 
or BNL energies [16] are of the order B ~ 10 17 — 10 19 G. In nature, we have another indication 
of sources of intense magnetic fields of astrophysical origin such as magnetars [17]. Magnetars are 
neutron stars which may have surface magnetic fields B ~ 10 15 G [18-20] discovered in the X-ray 
and 7-ray electromagnetic spectrum (for a review see [21]). They are identified with the anomalous 
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X-ray pulsars (AXP) and soft 7-ray repeaters [22]. 

Taking as reference the critical field, B^, at which the electron cyclotron energy is equal to 
the electron mass, B% = 4.414 x 10 13 G, we define B* = B/B^. It has been shown by several 
authors [23-26] that magnetic fields larger than B* = 10 5 will affect the EOS of compact stars. In 
particular, field-theoretical descriptions based on several parametrizations of the non-linear Walecka 
model (NLWM) [27] show an overall similar behavior. According to the scalar virial theorem [28] 
the interior magnetic field strength could be as large as B = 10 18 G so, in principle, this is the 
maximum field strength that is meaningful to consider. 

In 1959 the formal theory for treatment of low temperature (non-superfluid) fermion systems, 
known as normal Fermi Liquids, was developed by Landau [29, 30] to describe the behavior of 3 He 
below 100 mK. With this Fermi Liquid Theory (FLT) (see a recent reference [31]) the excited states 
in the system could be described as quasiparticles (qp) as long as these states have sufficiently long 
lifetimes. At low temperature, the small excitation energy (compared to the chemical potential) 
will assure this fact. In the context of the FLT, the so-called Landau parameters, can parametrize 
the interaction energy between a pair of qp in the medium. Previous works have attempted to 
partially study the behavior of Landau parameters for non-magnetized symmetric nuclear matter 
or neutron matter [32-34] or in magnetized matter under the presence of a magnetic field either in 
a non-relativistic formalism [35-37] or in magnetized matter without considering B field including 
exchange terms in a relativistic way [38]. 

In this work we will be interested in calculating the Landau Fermi Liquid parameters for an 
isospin asymmetric nuclear system in beta equilibrium and in charge neutrality under the effect 
of an intense magnetic field. The FLT used in this case must describe relativistically the more 
general condition of a magnetized non-pure isospin system to account for the fact that the intense 
magnetic field can modify isospin populations and partially align nucleon magnetic moments with 
respect to the case of vanishing magnetic field. In addition, the different dynamics of proton and 
neutron sectors under the presence of a magnetic field (including the existence of anomalous nucleon 
magnetic moments) will have effects in the Landau parameter computation showing discretized or 
continue features for protons and neutrons, respectively. In section II we introduce the relativistic 
lagrangian model used in this work and the generalized formalism of the FLT for charge neutral 
isospin asymmetric hadronic systems under the presence of a magnetic field. In section III we 
discuss the explicit form of the matricial structure of the coefficients describing the interaction of 
qp in the magnetized system through the Landau parameters. In section IV we analyze the obtained 
coefficients for either individual spin quantum numbers or total spin (singlet or triplet) for the qp 
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excitations for the electrically neutral system configurations calculated under beta equilibrium and 
the Landau parameter behavior in presence of a strong magnetic field and, finally, in section V, we 
summarize and draw some conclusions. 
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II. THE FORMALISM 

For the description of the EOS of neutron star matter, we employ a relativistic field-theoretical 
approach in which the baryons, neutrons (n) and protons (p), interact via the exchange of a — uj — p 
mesons in the presence of a uniform magnetic field B along the z-axis. The Lagrangian density for 
the TM1 parametrization [39] of the non-linear Walecka model (NLWM) reads [27] 

b=n,p m=a,Lu,p l=e 

The baryon (b=n, p), meson {m = <r, u and p) and lepton (/ = e) lagrangians are given by 

(c = h=i), 

C b = ^ b (i^ff* - -m b + g a a - g^^ - ^g p r^^ - ^N^a^F^A (2) 

-\f^F, u + \m 2 p p^ - \P^P, U , (3) 

Ci = ^>i{i ltl d^ -q ail A^ -mi)^ u (4) 

where *$> b , ^ are the baryon and lepton Dirac fields respectively. 

The nucleon isospin z-projection for the proton (neutron) is denoted by t^ p = 1 ( r^ n = — 1). 
The nucleon mass is m b {m b = m n = m p = 938 MeV), its charge is q b and the baryonic anomalous 
magnetic moments (AMM) are introduced via the coupling to the electromagnetic field tensor with 
a pu = | [jp,ju] and strength n b . In particular, K n = —1.91315 for the neutron and k p = 1.79285 
for the proton, mi and q\ are the mass and charge of the lepton. We will consider the simplest 
model where the leptonic sector is formed just by electrons (I = e), with no anomalous magnetic 
moment, providing charge neutrality in this astrophysical scenario. Despite there is a non-zero 
electron AMM its value [15] is tiny when compared to that in the hadronic sector and it was shown 
that this contribution is negligible for the magnetic fields of interest in astrophysics if properly 
introduced [40]. The mesonic and electromagnetic field strength tensors are given by their usual 
expressions: Q. pv = d p u u - d u uj p , P pu = d p p u - d u p pi and F pu = d^A u - d u A p . 

The electromagnetic field is assumed to be externally generated (and thus has no associated 
field equation) , and only frozen-field configurations will be considered in this work. To calculate the 
thermodynamic conditions in this charge neutral asymmetric nuclear system in beta equilibrium 
with an intense magnetic field, isoscalar and isovector current conservation must be imposed. 



Explicitly, the following conditions are fulfilled: i) electrical charge neutrality ii) conservation of 
baryonic charge iii) mesonic field equations selfconsistency. In addition, we will assume thorought 
this work that neutrinos scape freely and therefore there is no neutrino trapping. 

The field equations of motion are determined from the Euler-Lagrange equations arising from 
the lagrangian density in Eq.(l). Under the conditions of the present calculation, a relativistic 
mean field (RMF) approximation will be used so that the space-time varying fields are replaced by 
a homogeneous value, <f>(xu) — > <j>. In this way, the Dirac equation for a nucleon is given by, 

(«7m^ ~ Iblti^ - (m b - g a a) - g^l^ 

-\g P Tzbl^ - ^NWnvF^b = 0, (5) 

where the effective baryon mass is m* = m b — g a a. For leptons, 

(*7m^ - ~ m i) = °- ( 6 ) 

For meson fields we obtain, 



m\o + -reer 2 + ^Acr 3 = g a (n s p + n s n ) = g a n s , (7) 
m* w° + jtgi (oo° 2 - u 2 ) uj° = 9uJ (n« + <) = g u n B , (8) 
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+ gj^w (w° 2 - us 2 ) u = g w j B , (9) 



m2 pP° = \d P K ~ n n) = 7,9 P n 3 , (10) 

m pP = (11) 

where we use the notation as in the work of Matsui [34] and n B is the baryonic (vector) particle 
number density constructed as the sum of (vector) particle number density of protons (rip) and 
neutrons (raJQ, n B = Tip + n v n . 

The baryon current j B = jp + jn is also the sum of the proton (j p ) and neutron (j n ) currents. 
n s = rip + n s n is the scalar density constructed from that of protons (ni) and neutrons (n*). 
n 3 = n p ~ n n an d J3 = j p — j n are the isoscalar particle number density and isovector baryon 
current, respectively. When solving for equilibrium conditions in the nuclear system governed by 
Eqs.(5)-(11) we impose Jb = and j3 = 0. Then, we have for the a field, 

g a a = ^n s , (12) 
with m'^ = m 2 + \kcj + ^jAcj 2 . For the u° field we get 

gujOJ = -^n B , (13) 
to' 
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with m'l = ml + ^g^ 2 . 

When the Dirac equation for nucleons Eq.(5) is solved, a magnetic field B in the z-direction 
given by B = Bk is used. The energies for the quasi-protons and quasi-neutrons in the medium 
with spin z-projection, s, are given by the following expressions [23], 



m* 2 + 2vq p B — sji^KyB 



k 2 + [J m* 2 + k\- sfj, N n n B + g u u 



2 9pP 



->9 P P 



(14) 
(15) 



where v = n + k — sgn(qh)^ = 0, 1, 2, . . . enumerates the quantized Landau levels for protons with 
electric charge q p . The quantum number s is +1 for spin up and —1 for spin down quasipar- 
ticles. Due to the fact that the magnetic field is taken in the z-direction, it is useful to define 
three-momentum (k) components along parallel (k z ) and perpendicular (kj_) directions. Then, for 
neutrons, the surface of constant energy is an ellipsoid while, for protons, constant energy surfaces 
are formed by circumferences on nested cylinders with radius labeled by the Landau level, see next 
section. In this work we are mainly interested in hadronic properties, and electron dynamics will be 
such that at high B field strengths, they will mostly be in the low Landau levels. For completeness 
we write the expressions of the scalar and vector densities for protons and neutrons for both spin 
polarizations as follows [23] 
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where k F , k F are the Fermi momenta of protons and neutrons related to the proton and neutron 
Fermi energies, E P F and E F , respectively, by 

1 2 



K F,v,s 

h n2 
h F,s 



E p2 

F 



m* 2 + 2vq p B — sfj,NK p B 



E 



ii2 



m r . 



(17) 



and m n = m* — s[i]\rK n B = m s n . The summation in v in the above expressions terminates at 
u maxi the largest value of v for which the square of Fermi momentum of the charged particle is 



8 



still positive and corresponds to the closest integer from below denned by the ratio 



(E P F + sfj, N K p Bf 
%] B 



711 r 



, *2 



(18) 



For electrons fc|.„„ = \/ Ep 2 — (m 2 + 2vq e B) and f ma2: does not have the same value as that for 



protons. 



III. LANDAU FERMI LIQUID PARAMETERS 

We now calculate the Landau parameters using a generalized formulation of Landau Theory of 
Fermi Liquids [31], from the variations of the energy density of the system, £, given by [32] 

Iglml 2 I glmj 2 1 4 1 &g 4 2&« 2 2 

UJ UJ CJ UJ (jj 

lal o 1 9 2 



1 



1 At 



1 A 



*\4 



i? 2 

- ™*) 2 + ^r^-( m fe - m *f + T\~2^ mb ~ m *) 4 + 
2 5a 3! g l a 4! 5 £ 2 



(19) 



where we have defined n P = n p {ki,v, Sj) as the occupation number of the quasiprotons and n™ = 
n n (ki,Si) for the quasineutrons. For leptons the occupation number is n\ = n l {ki,u, Sj). Also we 
have defined the energies, 



<s = + ( \M 2 + toQpB ~ s^ n k p B 



(20) 
(21) 



and analogous for electrons, E ivs . We use generalized three-momenta depending on isospin 

Ki = ki-Vi, j=p,n (22) 



with V p = g w w + \g P p and V n = g^ — hgpP- The equation for the nucleon effective mass can be 
written as 



- {rip + re*), 



771 



(23) 



where 



S/ J,NK n B 

m* 2 + /C 2 ; , 



(24) 
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U ' P S ml s + Z N n p B' (25) 



where the following definition has been used 



m vs = Y m + 2q p uB - sk p ^nB = e p vs - SK p fiNB. (26) 
The vector current for quasiprotons is written as 



\ = V (27) 



and for quasineutrons, 



fc±i I -x sfi N n n B \ K n zi 



E n 

IS 



m* 2 + JCf, ) E i 



(28) 



so that ]b and J3 can be constructed. 

According to the FLT [31] the first variation of the energy density of the system, £ , with respect 
to the ocupation number for qp with isospin of jth-type, n^, defines the qp energy, e\. Let us notice 
that in reduced notation the index i means (i, v, s) for quasi-protons and (i, s) for quasi-neutrons: 

8£ = J24Sni (29) 

and 

^ = ^,. + ^ + J|" s ' (30) 

e ^ K + ^n B -^n 3 . (31) 

The single quasiparticle energy has, each, two explicit contributions, one due to the motion under 
the influence of a strong quantizing magnetic fiel and another due to the motion in a medium 
with mesonic self-interacting fields. To calculate the Landau parameters we use the standard 
approach [31] but generalizing to the case when there are external fields. These will allow to 
extract information on the interaction energies of quasiparticles of spin z-projection (s, s') and 
isospin and, if they are protons, different Landau levels {v, v') in the system. In a generalized 
system with a qp state the Landau parameters are calculated as the second derivative of the energy 
density of the system, £, with respect to the qp state with occupation number n\ gl , that is j-isospin, 
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spin s' and momentum ICj (if they are quasiprotons they include the additional quantum number 

l \vs' 



v, in the way nj ,). Using Eq. (29) it is defined, 



fL = Pr- (32) 

From the original formulation by Landau of the pure neutron system without considering spin 
degrees calculated in [29, 30] the above expression generalizes to a (2v max + 2) x (2v max + 2) 
matrix, fi, in isospin and spin space. In this way we have a characterization of the nuclear system 
when an external quantizing magnetic field is considered, 

(tpp f? n \ 
Jilss'uu' Jilss'u I (33) 
flip run J 

Jilss'v 1 •> ilss' ) 

The detailed calculation of the matrix elements is given in the appendix in section VI. 

In order to obtain the relativistic Landau parameters we must consider in the context of the 
FLT that the interactions of the effective quasiparticles in the system will take place close to the 
Fermi surfaces, since the lifetime of these excitations varies inversely with the departure of its 
energy, E, from the Fermi energy r ~ (e-Ef)' 1 ' 

Since the matrix elements in Eq.(33) have dimensions of energy divided by number density 
it seems convenient to define new dimensionless coefficients by multiplying them by the density 
of states at each Fermi level for quasiprotons, Np and quasineutrons, N%, at the Fermi surface 
with a given spin projection. In the case of protons the quantized level filling must be carefully 
considered and the definition of density of states at the Fermi level and in a given Landau level 
with polarization s is 



\ly\B p 

where \i = E V F is the proton Fermi energy. Summing over all possible levels we have, 



(34) 



where k p Fu s = y E p p 2 — (m^ s ) 2 '. For neutrons, 

'dn 1 ?\ E\ 



arcsm 



E F ) 2 



(36) 
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For the mixed states and due to the fact that the f lJ s , arise from two-body operators we can 
define the density, 

Nif = y/N?Nf. (37) 

Notice that this is not the only possible definition for the mixed density of states, however, the 
correct limit for isospin pure systems and B = limit can be recovered when used. Using this 
prescription, dimensionless coefficients F 1 ^, can be defined [29-31] so that: 

F*f=N!?f%. (38) 

The Fermi surfaces are selected between the energy surfaces given by Eq.(20) and Eq.(21) in the 
system at equilibrium. For protons these are cilindrical holes with k v z taking values ±k Fus at the 
Landau level defined by v for a given spin projection s as shown in Fig.l. For each spin projection, 
s = ±1 and a given v, the proton Fermi surface are two circles with kz = ikp . For neutrons, 
instead, the Fermi surface for given s is shown in Fig. 2 as an ellipsoid of maximum perpendicular 
momentum |fc|™'™ aa = kp ± s . This maximum Fermi momentum in the transverse direction is 

1^7 = k FXjS = ^{E F + s W B) 2 -m* 2 , (39) 

while in the z-direction the fc™ value attains a maximum value, 

\k\: : r X = k n Fs = y/Ef - (m-) 2 . (40) 

In regular FLT the relevant interaction takes place on the Fermi surfaces and then the Landau 
parameters depend on the density, ng, and the angle, 9, between qp three momenta being possible 
to perform their expansion in Legendre polynomials, Pi(cos(8)). Taking averages over angular 
dependence of the Landau parameters shows that the only terms remaining are the I = terms. 
Then, for a pure isospin system the relations // = l£i£^l±=± ; anc [ g ; = ^ s - s 2 ^ 3, ~ 3 hold and in this 
way values for Fq = NqJq and Go = N^go were the first to be historically obtained by Landau 
[29, 30] and can be obtained from indirect experimentally measured observables [34]. On general 
grounds we can define auxiliar combinations of Landau parameters, 
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and, from them, we define combinations with singlet (S = 0) or triplet (S = 1) qp states as, 



F (5=0) = F+ + FT. F {S=1) = G + + G~ 



(43) 



If the B — > limit is taken, we recover from the previous expression the usual Fq and Go in the 
normal FLT. 



A k 



c 



c 



p 



FIG. 1. Cilindrical shape of constant energy for quasiprotons. For spin projection s the Fermi surface are 
concentric circles with values k v z = k p F v and radius fhP s . Notice that as v grows k\ is smaller. 



In our work, in a similar way, and in order to evaluate the angular averages of coefficients 
over qp Fermi surfaces we must take into account the different dynamical behavior due to the qp 
electrical charge. For protons there are several Landau levels that can be populated over cilindrical 
holes with fef = ±kp^ s . The average of a function s, s', v, u') is performed as 

< f > S s'= ; ■ — — • (44) 

This definition gives the correct B — > limit as the one obtained with a regular FLT B = 
calculation. For neutrons the Fermi surface is an ellipsoid defined by Eq.(21) and the average of 
a given function g(k™,s) should be performed using the fact it presents axial symmetry. We will 
perform an integration over /c" resulting from the projection of the Fermi volume over the plane, 
that is a disk S± of radius \k\ n ^™ ax as given by Eq.(39). Then, in order to average we must replace 
the fej by its value on the Fermi surface 

k» = ±^Ef - (m?(k ± )) 2 , (45) 
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FIG. 2. Ellipsoidal shape of constant energy for quasincutrons. For spin projection s the Fermi surface has 
maximum values \k\f™ ax = k^ s and \k\^™ ax = k 7 * ± s 



with m^(kj_) = y fc" + m* — sfij^K n B, up to the maximum value, and integrate g(k™(k±), s) over 
the disk S± in the XY plane. The average of the function g finally reads as, 

IV. RESULTS 



We have used the thermodynamic conditions arising from the selfconsistent solution of the RMF 
set of equations Eqs.(7)-(11) solving beta equilibrium in a charge neutral nuclear homogeneous sys- 
tem under the influence of a strong quantizing magnetic field. In Fig. 3 we show proton population 
fractions for particles with magnetic moments polarized parallell, (solid line), and antiparallel, 
(Y~) (dashed line), to the magnetic field, B, for different baryonic densities ne/no = 0.1,0.5,2,4 
as a function of the logarithm (base 10) of the magnetic field strength. Due to the tiny value of 
the baryon magnetic moment, magnetic field strengths larger than Log B(G)> 16 are needed so 
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FIG. 3. (color online) Proton fraction for the up (red solid line) and down (green dashed line) spin as a 
function of the logarithm (base 10) of the magnetic field strength at baryonic densities ns/n = 0.1, 0.5, 2, 4. 

that there is a rapid increase (decrease) of the up (down) proton fraction. For a given magnetic 
field strength, the differences between fractions of protons polarized parallell and antiparallel to 
the magnetic field direction are larger for the smaller densities. This behavior is a consequence 
of the energy balance of the interaction of the proton orbit magnetic momentum and AMM with 
the magnetic field. Due to the positive charge and the fact that k p > 0, a quasiparticle energy 
is lowered by aligning spins with B, and the opposite for spin down. Let us remind that for den- 
sities below saturation density, n# = uq = 0.145 fm -3 , the spatially non- homogeneous systems 
are energetically prefered over the uniform ones with the onset of Pasta phases [11], however this 
density can be of interest to the study of low density neutron rich gas among the expected pasta 
structures. 

In Fig. 4 we show relative polarization, Aj = (nf — n~)/(nf + n~), where nf are the vector 
number densities of i-particle population, for neutrons (solid line) , protons (dashed line) and 
electrons (dotted line) as a function of the logarithm (base 10) of the magnetic field strength (in 
Gauss) at a density n# = 0.0145 fm~ 3 . At this low density and for fields Log (B(G)) > 16 
there is a complete alignment of the proton sector where they will all be in the n = Landau 
level, while the opposite charge of electrons force them to be in the antialigned state with n = 
for a slightly smaller field. Instead, the antipolarization of neutrons occurs due to the different 
interaction (k u < 0) of the neutron AMM with the magnetic field. In this sense, this low density 
case could somewhat illustrate the properties of the homogeneous neutron gas and it remains to be 
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FIG. 4. (color online) Proton (green dashed line), neutron (red full line) and electron (blue dotted line) 
relative polarization as a function of logarithm (base 10) of the magnetic field strength (in Gauss) for 
ns/na = 0.1. 

seen if the interaction of the B field and other effects [42] could prevent the formation of proposed 
exotic superfluid states in the interior of NS. 

In Fig. 5 we show relative polarization, A, for neutrons (red solid line) , protons (green dashed 
line) and electrons (blue dotted line) as a function of the logarithm (base 10) of the magnetic field 
strength (in Gauss) for n^/no = 0.1, 0.5, 2, 4. Increasing densities are depicted with decreasing line 
width. Hadrons and leptons show different behavior according to the associated sign and value 
of their magnetic moment. For low density at ub/uq = 0.1 protons show a total polarization for 
magnetic fields larger than B ~ 10 16 G while for electrons the polarization is slightly smaller and 
with opposite relative sign because of its negative electrical charge. From this figure it is seen 
that as density grows the magnetic field strength needed to polarize a given population fraction 
is bigger, due to competing effects in Fermi energies. The different behavior between protons and 
neutrons is due to the orbital magnetic moment contributing only for protons. In general, it is 
quite difficult to polarize neutrons, and A n remains always below A p (in absolute value). 

In Fig. 6 we show the averaged and normalized Landau parameters F?? given by Eq.(44) and 
Eq.(46) for the different combinations of isospin {i,j} = {n,p} as functions of the logarithm (base 
10) of the magnetic field strength (in Gauss) for: a) ub/tiq = 0.5 (left pannel) and b) n^/no = 2 
(right pannel). F pp , F pn , F np and F nn are depicted from top to bottom in each pannel. We can see 
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FIG. 5. (color online) Proton (dashed green), neutron (solid red) and electron (dotted blue) relative polariza- 
tion as a function of logarithm (base 10) of the magnetic field strength (in Gauss) for ns/no = 0.1,0.5,2,4. 

the different polarization components, (++),(H — ),( — h), ( ), in solid, long dashed, short dashed 

and dotted lines respectively. The peaks appearing in the dimensionless coefficients involving 
protons, that is, (F pp , F pn , F np ), are due to the fact the level density depends on the Landau level, 
n, and spin projection, s. Instead, for neutrons the behavior is smooth since Fermi surfaces are 
ellipsoids. We can see that the parameters F??' for the low density case on the left pannel are 
negative, signaling attractive interaction, while for the high density case they are always positive, 
or repulsive, for the magnetic field strengths considered in this work. Results for B field strengths 
with Logio(B(G)) > 18 should be taken with care and considered as an extrapolation. At low 
B the hopping behavior is smooth due to the large number of Landau levels populated and as 
B grows the number of levels decreases and the separation between them increases (see Fig. 1). 
For sufficiently large B, k p Fus= _ 1 = 0, and beyond that strengh the density of states, N~ and 
all components in the F pp except for the (++) are not defined. The same happens for F pn ,F np , 
since they involve the proton spin down component. For these components there is a symmetry 
by replacing simultaneously spin and isospin indexes, in this way, for instance, the F p ~ and F~ p 
are equal. For strengths of B ~ 4 x 10 18 G all the neutrons are polarized down and the only non 
vanishing components are those involving this fraction. 

At low densities and for components with mixed spin the quasiparticle states are more bound 
than same spin components. Protons tend to polarize up, as B grows, filling the low Landau levels 
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FIG. 6. (color online) Averaged and normalized Landau parameters F^ s as functions of the logarithm of the 
magnetic field strength for s,s' = ±1 at densities rts/rio = 0.5 (left panncl) and ns/rtg = 2 (right pannel). 



and the qp interaction becomes more attractive, however for neutrons this behavior reverses since, 
the larger B, the smaller the neutron fraction populates the system and the more repulsive the 
medium reacts to the formation of a qp state. Then, this causes quasineutron excitations to be less 
bound on an absolute scale. In the high density case (right pannel), the spikes signal the hopping 
of the Landau levels, and the clear more repulsive interaction than shown on the low density case 
(left pannel) where the magnetization is stronger. The mixed states involving protons are less 
repulsive than the F nn in the mostly neutron populated system. 

In Fig. 7 we show the F$ normalized Landau coefficients according to the definitions in Eq.(43) 
with a definite value of total spin, S, for all combinations of isospin {i, j} = {n,p} as a function of 
the logarithm (base 10) of the magnetic field strength (in Gauss) for a) ub/uq = 0.5 (left pannel) 
and b) ns/no = 2 (right pannel). F pp , F pn , F np and F nn are depicted from top to bottom in 
each pannel. The singlet (solid line) and triplet (dashed line) state parameters present a different 
behavior due to their construction. As seen in the left pannel for the low density case, due to 

cancellations of negative components (++,H — , — h, ), the triplet states have small binding. 

Notice that for the F pp the well-defined values are those where k p Fus is real [y < v max ). For singlet 
states the energy is lowered with respect to triplet states, for different isospin components (pn, np). 
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Log 10 [B(G)] Log 10 [B(G)] 



FIG. 7. (color online) Averaged and normalized Landau parameters as functions of the logarithm of the 
magnetic field strength for S = 0, 1 for densities ns/no = 0.5 (left pannel) and ub/tiq = 2 (right pannel). 

At higher densities, ub/tiq = 2 (right pannel), the different components F?? are all positive, 
and, therefore, the singlet and triplet binding behave in a similar way to the low density case, 
however, in absolute value the singlet case will have repulsive character. This analysis of Landau 
parameters could be of potential interest to obtain transport coefficients from the components of 
their expansion in cylindrical harmonics, due to the axial symmetry introduced by the existence of 
a strong quantizing magnetic field. Prom this expansion information of low temperature nuclear 
system observables [41] such as, for example, effective mass, spin diffusivity, viscosity, etc could be 
accessed. Additional contributions from exchange terms to the calculated parameters should also 
be considered. 

V. SUMMARY 

In the present paper we have derived for the first time the Landau parameters for charge 
neutral homogeneous asymmetric nuclear matter under very strong magnetic fields. In particular, 
we have considered cold stellar matter in /3-equilibrium in the presence of fields with strengths 
Logio (B(G))< 18 as allowed by the scalar virial theorem. We have used a RMF model to describe 
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the nuclear system so that our calculation may be extended to high densities. However, no exotic 
matter such as hyperons or kaon condensates has been considered. In our work, we have only 
considered matter densities from subsaturation density, ub/uq = 0.1, to densities before a possible 
quark transition. Present results at subsaturation densities should be taken with care, as a first 
approach to a realistic description of the effect of strong magnetic fields on the cold stellar matter. 
Also we have not considered the existence of superfluid neutrons at subsaturation densities. 

We find that the nuclear system needs at least magnetic field strengths larger than B ~ 10 16 G 
to show some noticeable magnetization in excitation spectrum. Protons tend to polarize aligning 
spins with B field direction, while neutrons and electrons do the opposite. Landau coefficients F*? 
show the interaction energy of qp excitations with isospins i and j and spins s and s', respectively, 
and are calculated from the energy density £ of the system in charge neutrality as a two body 
operator obtained from its second derivative with respect to the occupation numbers of the qp 
excitation. In the excitations involving quasiprotons the discrete Landau levels play an important 
role quantizing the momentum states in the direction perpendicular to the field. In this way there 
is a discrete feature associated in the Landau coefficients. For small B the number of allowed levels 
is large and is hardly visible on the plots, however for higher B the number of levels decreases 
and they are more separated. For neutrons instead, the lack of quantized Landau levels implies a 
smooth behavior of the parameters as functions of magnetic field strength. 

At low densities the effect of the magnetic field is small on the individual polarization since 
the polarization induced is weak (for fields with Logio(-B(G)) < 17) and gives an attractive overall 
effect with respect to unmagnetized systems. The dominant component drives the tendency of 
the parameter, in this way the proton sector is more bound and so are the pp,pn,np components, 
while for neutrons they are less bound since the number of neutrons decreases when B increases. 
In this case, for situations of complete polarization the Landau parameters get largely negative. 
In the high density cases the interaction is increased in the way to unbind the qp excitations. The 
medium effects are such that the interaction is repulsive and the magnetic field tends to decrease 
this effect in the less populated fractions while increasing in the more populated ones. 

For the singlet qp configuration the medium effects play a role of binding at low densities while 
at high densities the opposite occurs. The B field tends to reverse this tendency (if not too strong) 
and can be considered a competing effect. As for the triplet the qp spin configurations the in- 
medium effects play no role, while the B field brings extra binding in the pp or nn channel at low 
densities, while the opposite occurs at high densities. The effect of the B field is reversed in the 
pn or np channels, binding (unbinding) is present at high (low) densities. 
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The calculation we have performed could have potential interest to size the relevance of the 
inclusion of a magnetic field in astrophysical scenarios with realistic description of nuclear systems 
in beta equilibrium by relativistic mean field theoretical models. We have obtained singlet and 
triplet interaction energies of qp states and further work could relate these, in a consistent fashion, 
to relativistic and (possibly) superfluid hadronic components in this type of systems. 
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VI. APPENDIX 

In this appendix we provide the explicit form of the quasiparticle interaction matrix elements 
depending on spin and isospin and the mass and current derivatives appearing in the calculation 
of the Landau parameters in section III. 

A. Interaction matrix elements 



The matrix elements appearing in /x, in Eq.(33) in section III read as: 
i) for quasiproton-quasiproton interaction, 

del. . 



fpp 

J ilss'vv' 
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(47) 



ii) for quasineutron-quasineutron interaction, 
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(48) 



iii) for quasiproton-quasineutron interaction, 
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A. 
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(49) 



iii) for quasineutron-quasiproton interaction, 
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The indexes in labels in expressions Eqs.(47)-(50) have been accordingly written with six numbers 
(i, is, s), (I, v' , s') for quasiproton interaction, and four numbers (i, s), (I, s') for quasineutron inter- 
action. For isospin mixing terms (quasiproton (quasineutron)-quasineutron (quasiproton)) there 
are 5 numbers (i,v,s), (l,s f ) in the label. We have used auxiliar definitions in the combination of 
coupling constants, 

,2 n 2 



A± = \± 

m„f 



4 mi 



and the proton effective energy 



(51) 



(52) 



with an auxiliar mass m v vs = ^m* 2 + 2q p vB — sk p [inB = e£ s — sKpfi^B. 

For neutrons we define an auxiliar energy (which can be interpreted as a neutron momentum 
dependent mass): 



e±i = \l m *n + k 2 



and the effective energy for neutrons: 



(53) 



(54) 



We also define, 



D = (1 + A + b n ) (1 + A + b p ) - A 2 _ b p b n 



where the coefficients b p and b n are: 
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and k Fs 



El 2 



m 2 . The explicit form for the integrals I±, I2 and I3 is given by 
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The integral I2 can be written out into several pieces l<i = 1% + l\ + I$- 
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with q = y / (Eg) 2 — (sK n /j,]\fB) 2 . The integral ^3 can be written as, 
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We write the integral as several pieces 1$ = If + I3 + If. For If we have = J|. For l|, 
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B. Effective mass derivatives 

The mass derivative appearing in eqs. (47-50) — = 9 (™^~9°< T ) can De obtained using the 

Is' ' v 1 Is' v' 

equation of motion for the scalar a field in Eq. (7). Then we can write 
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In the limit where sums are converted into integrals it is found that, 
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with the integrals I\ and I2 having the form given in Eq.(57) and Eq.(59). 



C. Current derivatives 

The baryonic currents Eqs. (27)- (28) in the asymmetric nuclear system can be written in terms 
of the components in directions parallel k and perpendicular h± to the external magnetic field. 
Using explicitly the expression for the effective momentum K,\ = fcj — V^, j = p, n we have, for 
protons 

j P = jpzk, (74) 

with 
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and for neutrons 
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Using the conventions as defined in Matsui's work [34] j# = j p +j n an d j3 = j p — j n the derivatives 
of the baryonic isovector current components are obtained assuming that the macroscopic currents 
in equilibrium will vanish, i. e., j# = and j3 = 0. In this way the derivative of the baryonic 
vector current (z and _L-components) with respect to the proton density nf are: 
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with integrals J\ and J 2 that can be decomposed as J\ 
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For the derivatives of je with respect to the neutron density nf we have 
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with the integral J3 which can be written out as J3 

<9j 3 | 2 r/c? 



h. If we now calculate the derivatives of 



the J3 current with respect to nf 
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and with respect to neutron density nf , 
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To solve the coupled set of Eqs.(77)-(84) for the derivatives of jb and j3 the following formula 
can be used in a straightforward way. For the system, 
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